Noiseless nonreciprocity in a parametric active device 
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Nonreciprocal devices such as circulators and isolators belong to an impor- 
tant class of microwave components employed in applications like the mea- 
surement of mesoscopic circuits at cryogenic temperatures [1, 2, 3, 4]. The 
measurement protocols usually involve an amplification chain which relies on 
circulators to separate input and output channels and to suppress backaction 
from different stages on the sample under test. In these devices the usual recip- 
rocal symmetry of circuits is broken by the phenomenon of Faraday rotation 
based on magnetic materials and fields [5]. However, magnets are averse to 
on-chip integration, and magnetic fields are deleterious to delicate supercon- 
ducting devices [6, 7, 8]. Here we present a new proposal combining two stages 
of parametric modulation emulating the action of a circulator. It is devoid of 
magnetic components and suitable for on-chip integration. As the design is 
free of any dissipative elements and based on reversible operation, the device 
operates noiselessly, giving it an important advantage over other nonrecipro- 
cal active devices for quantum information processing applications. 
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Reciprocity is one of the fundamental symmetries frequently encountered in electrical cir- 
cuits. It is equivalent to the more familiar notion of the principle of reversibility in optics which 
states that any experiment is symmetric under an exchange of source and image [9]. Reci- 
procity can, however, be violated, for example, by the magneto-optic effect of Faraday rotation 
[5] which leads to rotation of the polarization vector of light resulting from different propaga- 
tion velocities of left- and right-circularly polarized waves in the presence of an applied mag- 
netic field B parallel to the direction of propagation (Fig. 1). The nonreciprocal phenomenon 
of Faraday rotation should be contrasted with the superficially similar, though reciprocal, ef- 
fect of optical activity where the polarization vector of light is rotated on passage through a 
non-centrosymmetric (chiral) medium. This change in the sense of rotation of polarization for 
counterpropagating waves in a Faraday medium (as seen by the observer receiving the light) 
has led some physicists to refer to this effect as a form of time-reversal-symmetry-breaking, a 
use of words that we prefer to avoid here [10, 11]. 

The phenomenon of nonreciprocity has propelled numerous theoretical investigations [[12] 
and refs. therein]; furthermore it offers immediate practical applications. Recent progress in 
solid state superconducting qubits, that provide some of the most promising architectures for 
scalable quantum computers [13, 14], has generated a huge incentive to integrate the compo- 
nents required for qubit operations and readout on-chip for incorporation in future quantum 
mechanical processors. A large variety of qubit readout protocols involve microwave reflec- 
tion based measurements and rely on nonreciprocal devices like circulators (or isolators) for 
separation of input and output channels [1,2]. These devices also play a strategic role in mea- 
surements based on low-noise microwave parametric amplifiers which, with the exception of 
designs based on the current-biased dc SQUID (Superconducting Quantum Interference De- 
vice), are also operated in a reflection mode with both the input and output signals collected 
on the same spatial channel [15, 16]. However, circulators (and isolators) routinely use bulk 
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components made from ferrites to achieve nonreciprocal phase shifts (Fig. lc) through Faraday 
rotation, making them unamenable to chip fabrication. Moreover, to bias the magnetic field 
in the ferrite, most of these devices use a permanent magnet which may channel flux into the 
superconducting device under test. 

In this letter, we present the full analysis of a model for a four-port circulator based on 
parametric active devices with no magnetic components. In active devices the energy source 
- provided by the pump - acts as the external "bias" field and sets the reference phase for the 
system, in analogy with the role played by the magnetic field in a Faraday medium. We exploit 
this effect in a cascade of active devices with pump phases at each stage tuned appropriately to 
obtain nonreciprocal transmission. 

The main building block of our design is a reversible IQ (in-phase/quadrature) modulator 
capable of performing noiseless frequency up- and down-conversion. A convenient analytical 
model capturing the fundamental properties of the device is shown in Fig. 2a. The device com- 
prises two low frequency LC resonators (addressed by two semi -infinite transmission lines A 
and B) coupled to a high frequency resonator (addressed by the transmission line C) through 
time- varying couplings M x , M 2 that emulate the role of the pump drive in active nonlinear 
devices and transfer energy from the tone at lu c to the signal modes propagating on the trans- 
mission lines. It operates in a manner analogous to the IQ modulation schemes routinely used 
in radiofrequency (RF) communication systems and microwave pulse engineering (hence the 
name) and converts two orthogonal spatial modes travelling on two distinct spatial ports (A, B) 
at same frequency (here uoq) into two orthogonal temporal modes travelling on the same spa- 
tial line (C) at different frequencies (o>+,u;_). In view of the reversible frequency conversion 
performed by this device (Fig. 2a), we will henceforth refer to it as the up/down-converter 
(UDC). In practice, such a device can be implemented on-chip using a ring modulator based on 
Josephson junctions, along the lines of the recently demonstrated experiment with Josephson 
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parametric converter [17, 18]. 

The complete design for the active circulator (Fig. 3a) consists of a UDC functioning as a 
frequency up-converter, a phase- shifter and a second UDC functioning as a frequency down- 
converter. 

A concise representation of the dynamics at each of the three stages in the cascade is pro- 
vided by the scattering matrix S which relates the outgoing wave amplitudes to the incoming 
wave amplitudes as seen from the ports of a network. We start by deriving the scattering matrix 
of the UDC stage. This is done by evaluating the impedance matrix Z of the UDC, as seen from 
its ports, and using the identity [5] 



S = (Z + Zo)' 1 x(Z- Z ) 



(1) 



where 

Z = diag(Z j4 , Z B , Z c , Z c ), (2) 

with Z A = Z B and Z c denoting the characteristic impedances of the semi-infinite transmis- 
sion lines serving as low and high frequency ports respectively. We obtain (see supplementary 
information, Fig. SI) 
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Here a and b denote the (reduced) amplitudes or the annihilation operators for the waves travel- 
ling on left and right transmission lines respectively (see supplementary information for details). 
These satisfy bosonic commutation relations of the form [19] 



[a h a\] = 5(uji -Uj). 



(4) 



In writing Eqs. (3) and (4), we have set a = a[u ],a + = a[a; + ],a_ = a[a>_] (see Fig. 2b). 
Similarly, the reflection coefficients at various ports are denoted by r ,r + and r_. The cross 
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reflection between the low frequency signal ports is denoted by q . The transmission coefficients 
are written as t (transmission without conjugation) and s (transmission with conjugation) with 
subscripts (u, d) indicating the up-conversion and down-conversion respectively. It is useful 
to note that the phase of the carrier, denoted by 0, affects only the transmitted amplitudes 
and rotates the two sidebands in opposite directions as can be seen from the corresponding 
scattering coefficients s and t in Eq. (3). The invariance of reflection amplitudes to the phase of 
the coupling will be important in understanding total reflections of the cascade, as we describe 
later. 

Further, we note that the matrix obtained in Eq. (3) is non-unitary, that is S^S ^ 1, which 
implies nonconservation of photon number as is natural for an active device. The matrix re- 
covers its unitary form as we turn off the couplings M 1? M 2 responsible for energy transfer 
between the pump and the signal modes. The full 8x8 matrix [supplementary information, Eq. 
(SI 8)], describing the device operation for all modes and their respective conjugates, fulfils the 
fundamental requirement of symplecticity\ll]. 

We can similarly describe the action of the frequency-independent phase shifting (PS) stage 
using a scattering matrix of the form 
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For each stage of the cascade, we now go from the scattering matrix representation to the 
transfer matrix representation [5], 
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since it is straightforward to calculate the total transfer matrix of the device by multiplying the 
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respective transfer matrices of different stages [20], 

^total = T DCr x T PS x T UCl ■ (7) 

Here the subscripts L, R index the left hand upconversion (UC) and right hand downconversion 
(DC) stage (Fig. 3a). The scattering matrix of the whole device is then obtained from T tota i 
using the inverse of the transformation in Eq. (6) (see supplementary information). We also 
note that 

T DC = F-'xT^xF 

= FxT^xF, (F- X = F) (8) 

where F — a x ® h, (&X is the 2D pauli spin matrix and I 2 is the 2D unity matrix). This matrix 
F is required to flip the indices, thus maintaining consistency in labelling the 'in' and 'out' 
amplitudes along a given direction of propagation. 

In our analysis we consider the operation at resonance, that is, when the input signal fre- 
quency coincides with the band center of the input resonators. Setting the phase of the pump at 
the first UDC stage 4>l = for calculational simplicity, we observe a transmission resonance 
for 9 = ±7r/2 (phase rotation by the PS stage), <f) R = ir/4 (phase of the pump at the second 
UDC stage), 5± = l/y/2 (detuning of the sidebands from the carrier in units of linewidth i.e. 
half width at full maximum of the resonance lineshape) of the high frequency resonator), and 
«l = = M /^Jl AjB L c = 2~ 3 / 4 (strength of the parametric coupling). For this choice of 
parameters, we obtain the scattering matrix of the complete device as 
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This is the matrix of a perfect four-port circulator. The analogy between a conventional circula- 
tor and the active circulator design proposed in this paper is made apparent from the respective 
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wave propagation diagrams in Figs, lc and 3b (see supplementary information for details on 
the calculation of coefficients on different arms in 3b). Nonetheless there are important differ- 
ences between the two designs despite the identity of the final S matrix. The coefficients on the 
forward (green) and backward (red) propagating arms of the active circulator design (Fig. 3b) 
involve deamplification followed by amplification, unlike the passive splitters (90 or 180 degree 
hybrids) employed in Faraday rotation schemes. This can be observed by squaring the ampli- 
tudes on each of the two arms originating from (or terminating into) a port and calculating the 
net power output, for each isolated UDC stage. It is straightforward to observe that, unlike the 
case of Fig. lc, they do not add up to unity. Nonetheless, the overall transmission is unity due to 
an exact cancellation of the reduction and gain in amplitudes. The wave propagation diagrams in 
Fig. 3b reveal another important difference of this design from that of a conventional circulator. 
The non-reciprocal action of the active circulator is not based upon any non-reciprocal phase 
shifters; instead it relies on the active stages used for frequency up- and down-conversion. The 
phase matching condition in the forward direction is met by tuning the phase of the coupling 
at the input and output UDC stages. In the reverse direction the phase mismatch leads to unity 
transmission in the spatially orthogonal port instead, leading to complete isolation between the 
incident signal port and its corresponding output port. 

Fig. 3c shows a convenient method to visualize this circulator action geometrically by 
mapping the device dynamics at different stages using a modulation ellipse. This approach is 
inspired by the polarization ellipse used to represent of state of polarization of an electromag- 
netic wave (linear, circular or elliptical), which involves recording the trajectory traced out by 
the tip of the polarization vector of light (defined by the instantaneous direction of the electric 
field vector E) in a plane perpendicular to the direction of propagation. Equivalently, a two- 
dimensional representation of the components E x and E Y of the electric field in the complex 
plane can be used to obtain a geometric description of the polarization of the light wave. In 
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the case of a modulation ellipse representation of the dynamics of the proposed device, we ex- 
tend this idea to map two distinct orthogonal modes (x,y) at each stage of the device [spatial: 
(x = a' ,y = Oq) or temporal: (x = a + ,y = a_)] as an ellipse in the plane defined by the 
coordinates / = Re[x + y\ , Q = Im[x — y*\ . This exercise shows that the final ellipses obtained 
at the output in case of forward and backward propagation through the device are rotated by 
90 degrees with respect to each other. This indicates that in the case of reverse propagation 
the orthogonal spatial port, relative to the forward propagation, receives the transmitted energy 
leading to a circulator action (see supplementary information and Fig. S2 for more details). 

Furthermore, as seen from Eq. (3), the reflection coefficients at the UDC stages are non- 
zero for all modes. However, for the whole cascade, the total reflection is identically zero at 
every port [su = for all % in Eq. (9)]. This remarkable cancellation of total reflections for 
the cascade can be understood in analogy with a Fabry-Perot resonance where a cavity flanked 
by two identical reflecting mirrors displays unity transmission at resonance. The total phase 
shift between the active "mirrors" in our device: (7r/2) a+ — (— 7r/2) a t = n, is akin to the 
resonance condition when a half-wavelength of the incident radiation equals the length of the 
Fabry-Perot cavity. Also the reflections at the two UDC stages are identical [as the reflection 
coefficients are independent of the phase angle 0, cf. Eq. (3)], fulfilling the second condition 
for the transmission resonance and net cancellation of reflections [21]. 

We show the dependence of circulator action on different parameters in the device in Fig. 
4. Since the isolation achieved is robust to reasonable deviations of parameters from their ideal 
values ((f) L = 0, 4>r = vr/4, 9 = 7r/2, 5± = l/y/2, cti = (*r = 2~ 3//4 ), the active circulator 
design holds promise for use in practical circuits. Another interesting feature of this device 
is the reversal of transmission characteristics with phase of the pumps (4>l,r —4>l,r) (Fig. 
4c). In the classic circulators based on passive Faraday rotation, this can be accomplished by 
changing the polarity of the magnetic bias field. Thus the clock ("pump phase") in an active 
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device indeed plays a role equivalent to the magnetic field in a Faraday medium. 

In conclusion, we have described a scheme for achieving nonreciprocal wave propagation 
using a protocol involving up-conversion followed by down-conversion mediated by an appro- 
priate phase shift [22]. The proposed design performs noiselessly as it consists of purely disper- 
sive components with no dissipation, making it attractive for quantum information applications 
using superconducting circuits [23]. Besides microwave applications, the architecture of the 
protocol described in this paper can be adapted to optical frequencies, where it can complement 
the recently proposed designs of nonreciprocal light propagation based on dynamical modula- 
tion of the refractive index of photonic structures [24] and the use of a surface waveguide on 
photonic crystals [25]. In addition to the practical applications outlined above, the treatment de- 
scribed in this letter may also give theoretical insights into the inherently directional dynamics 
of devices like the dc SQUID [26], when additional active stages are included in the chain. This 
can be useful in tackling unanswered questions pertaining to the quantum noise of dc SQUID 
amplifiers [27]. 
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Figure 1 : Faraday rotation and circulator action, a, Faraday rotation for a wave travelling 
from left to right in a Faraday-active medium, followed by a reflection back into the medium 
leading to a reversal of the direction of propagation. The rotation of the light polarization 
is fixed to a rotation-like property of the medium (shown by the arrows), set by an external 
magnetic field oriented along the propagation axis. The sense of light rotation as seen with 
respect to the direction of propagation remains the same, leading to the doubling of the rotation 
angle on reversing the ray through the medium, b, Rotation of the polarization vector of light 
on passage through an optically active medium, on the other hand, cancels out on reversing 
the direction of propagation. This occurs because optical rotation depends on the chirality of 
the medium (represented as a helix) which also reverses with the direction of propagation, c, 
Representation and schematic design of a conventional four-port circulator. The device consists 
of two 90 degree hybrids (equivalent to optical beam splitters) separated by a non-reciprocal 
phase shifter based on Faraday rotation. Solid black arrows indicate an amplitude split with no 
phase change while, open arrows indicate an amplitude split with a 90 degree phase change. 
The non-reciprocal phase shift is effective only for the propagation direction indicated by the 
arrow on the phase shifter box. 
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Figure 2: Description of an active reversible (information-conserving) IQ modulator per- 
forming frequency up- and down-conversion (UDC). a, Circuit schematic of the UDC 
containing only dispersive components. The two low frequency series LC resonators (with 
La = Lb and Ca = Cb), are fed by two input semi-infinite transmission lines, A and B, 
and parametrically coupled to a third high frequency series LC resonator leading to an out- 
put line C. The parametric coupling is achieved by varying the mutual inductances Mi and 
M 2 between the left and right resonators at the carrier frequency co c which, for optimal fre- 
quency conversion, is set at the band center of the right resonator. When operated from left 
to right, the circuit performs the modulation of low frequency signals of frequency ojq travel- 
ling on ports A and B to generate sidebands at co c ± co travelling on the high frequency line 
C. It performs the inverse operation of demodulation when operated in reverse from right to 
left, b, Spectral density /response landscape for different spatial channels (or ports) of the UDC 
circuit in a as a function of frequency. The dotted lines represent the couplings between dif- 
ferent ports. The solid and the dashed arrows represent different frequencies and respective 
conjugates. The resonance lineshapes of the two spatially distinct ports A and B are centered at 
oo a = ojb = V \J La,bCa,b- Here we show the case when the incoming signal at co is resonant 
with the center frequency (uj = uja,b)- The two sidebands generated by the UDC on channel 
C are detuned from the carrier u c by equal amounts. 
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Figure 3: Description of the active circulator, a, Circuit schematic of the active circulator 
design: the first UDC stage acts as a frequency up-converter (UC) (also indicated by a gradation 
in the color of the relevant box) with a parametric coupling modulated at the carrier frequency 
ui c = cu + — cu = <jj- + co and a phase 0t = 0. This is followed by a phase shifter (PS) that 
phase shifts both the sidebands by tt/2, in opposite directions. They are then demodulated by 
the final UDC stage acting as a frequency down-converter (DC), with the carrier phase 4>r = 
7r/4. b, Forward (green) and backward (red) propagation diagrams calculated using transfer 
matrix method for a with appropriate choice of detuning (S± = 1/ v2) and coupling strengths 
(oc L = a R = 2" 3 / 4 ) for maximum isolation, c, Representation of the device operation using 
modulation ellipses at each stage in the cascade. The top panel shows the forward propagation 
for the case when two distinct signals enter ports 1 and 2 respectively (a = 1, a = 0.5e" ? ; 
r\ = 7r/24) while the bottom panel shows the backward propagation dynamics when the same 
signals enter ports 3 and 4 respectively (b = 1, b = 0.5e ir? ). The relative phase and amplitudes 
are chosen to represent the most general case of two input signals which differ in both amplitude 
and phase. The relative phase difference between the two signals is encoded as the tilt of the 
modulation ellipse in the IQ plane while their average phase is represented as the color along the 
perimeter of the ellipse with yellow indicating zero phase (also see supplementary information 
for more details). 
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Figure 4: Variation of the difference between forward and backward transmission coeffi- 
cients (|s3i| — |si 3 |). Asymmetry in transmission, calculated for coupling angles <fi L = and 
4> R = re/ 4, as a function of a, strength of the coupling a and phase rotation 9 performed by the 
second phase shifting stage, b, detuning 5± of the sidebands from the carrier and phase rotation 
9, and c, detuning S± and the phase of the pump at second UDC stage 4>r. The points of maxima 
correspond to the ideal values reported in the text. The plot in b, also shows the periodicity of 
the response of the device as a function of 9. In c, the variation with respect to the pump phase 
shows the reversal of transmission characteristics with 0# — >■ —<f)R. As in b, the response is 
periodic in 4>r with a period equal to it. It can be seen that the design continues to work for 
moderate deviations from the preferred phase angle 9 = n/2, coupling otL,R = 2 _3/ ' 4 , detuning 
5± = 1/ a/2 and pump phase 4> R ■ = tt/A (values indicated with dashed arrows along the axis). 
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Methods 

-a 
a 

Derivation of scattering matrix. In this section, we present the details of the derivation of 
the scattering matrix S for the U DC stage. A scattering matrix has the general form 

a[u] out = S[u,oj']a in [Lu'], (SI) 

where oW 011 * is a column vector formed by the reduced mode amplitudes (in terms of photon 
number at the relevant frequencies) of the network. This description extends to all the spatial 
and temporal modes of the network. A convenient and simple way of deriving the S matrix 
involves calculating the impedance matrix Z using usual circuit theory, and then using the 

a 

identity [1] 

s = (z + z y\(z-z ) (S2) 

to obtain the S matrix. Here Z denotes the characteristic impedance of the transmission lines 
acting as channels for propagation of the incoming and outgoing signals at various ports. It is 

1 




RF output 



Figure SI: Comparison of two IQ modulation schemes: a Detailed circuit schematic of the 
UDC performing reversible IQ modulation. The fluxes across different inductances denoted 
as $j and the respective mode currents denoted as Ij are shown for each series LC resonator. 
The mutual inductances Mi , M 2 lead to off-diagonal coupling (mixing) between the current of 
the left resonators and the flux across the inductance of the right resonator and vice versa, b, 
Block diagram of a non-reversible circuit performing the operation described in a. The circuit 
employs dissipative components like the mixers (represented by the cross in circle symbol). 
The CW tone from the generator imposing a carrier frequency tu c is split into two copies, one 
phase shifted by n/2 with respect to the other. The carrier is modulated by the signal at ujq and 
encoded on two separate channels: the 'in phase" (or /) component and the "quadrature" (or Q) 
component. The / and Q channels are then combined to propagate onto a single spatial channel. 

of the form 

Z = diag(Z j4 , Z B , Z c , Z c ), (S3) 

where Za = Zb and Zq are characteristic impedances of the semi-infinite transmission lines 
addressing the low and high frequency ports respectively. 

For the parametrically coupled series LC oscillators forming the UDC stage, Fig. SI a, we 
obtain the total impedance matrix Z by adding the inductive (Z\ n &) and capacitive (Z cap ) contri- 
butions respectively. The inductance matrix L defines the constitutive relationship between the 
currents and fluxes for different inductances of the circuit: 
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In writing the above matrix, we ignore the fluxes at higher harmonics of the modes at cu and 
cu±. The inductive contribution to the impedance is then calculated by using the identity 

(S6) 

V/-/ 
V[u>] = Z iQd I[cu] 

where the subscripts denote the relevant frequency modes. It is straightforward to define the 
capacitance matrix for the circuit in the same manner: 
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The extra negative sign in s 44 of Eq. (S7) accounts for the generation of the conjugate wave 
amplitude a^[u;_] as a result of mixing of carrier u c and signal uj . Hence, on taking the Fourier 
transform of the current and voltage vectors, an extra negative sign appears for the correspond- 
ing coefficient in Z cap . The total Z matrix can then be written as 
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which gives 
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with Ti denoting the linewidth of the i th resonator. On using Eqs. (S9) and (S3) in Eq. (S2), we 
obtain 
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where a^ out denotes the shorthand notation for a[ui] u ^ out [cf. Eq. (SI)]. The symbols a, b 
denote the wave amplitudes for left and right ports respectively. The detailed expressions for 
scattering coefficients are listed below: 
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■<c / \(5 + i)(5±+i) + 2a 2 ) ' 
Equations (S12)-(S17) show that the effective coupling strength a plays the role of the paramet- 
ric drive ("pump") in the UDC. In the limit a = 0, transmission coefficients U and Sj become 
identically zero while the reflection coefficients r\ reduce to those for three independent series 
LCR circuits with resonance frequencies lu and uo c respectively. 
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We note that the matrix in Eq. (Sll) is one of the block diagonals of the full 8x8 scat- 
tering matrix that describes the interaction of both the positive and negative frequency wave 



amplitudes for all participating modes: 
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This is especially important in the case of an active network where a conjugation operation is 
possible i.e. Oj h-> a\ = a[—Ui\. The full scattering matrix S satisfies the following general 
properties: 

1. det(S)=l (as det (S) =1) 

2- X^=l \ s ij\ 2 — 1 

3. S T JS = J, 

where J represents a symplectic structure defined on the 2N x 2N phase space (N = number 
of degrees of freedom), 

J = ia Y ®I N . (S19) 

The last condition of symplecticity follows from the fact that a transformation of the modes as 
performed by the scattering matrix needs to be a canonical transformation. This requirement 
translates into the condition for preservation of phase space volume (information) or the number 
of participating modes (degrees of freedom) in the system. 



Derivation of transfer matrix. In this section, we derive the transfer matrix for the UDC 
stage. This description is equivalent to the usual ABCD matrix of the circuit theory defined in 
terms of the voltages and currents for a two-port network [1], 
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As there exists a straightforward mapping between the reduced wave amplitudes a, introduced 
earlier, and the currents and voltages at the ports [2], 

Vi + Zol, 
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we can easily adapt the concept to the above choice of variables. 

For the cascaded chain, we first evaluate the transfer matrix of each stage and then multiply 
them to obtain the total transfer matrix of the cascade. For the UDC stage performing up- 
conversion, we can find the relevant transfer matrix by solving Eq. (SI 1) to obtain sideband 
amplitudes (a+, at) in terms of low frequency amplitudes (a' Q , %) 
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The reversal of 'in' and 'out' in the column vectors on left and right hand sides of Eq. (S23) is 
required to maintain a consistent sense of propagation through the device as the output of the 
(N — l) th stage acts as the input for the N th stage in the chain. On doing the above transforma- 
tion, we obtain 
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The subscripts (+,-) refer to the resultant sideband generated at the output (u + or cu*_) while 
LR(RL) indicates the relevant direction of propagation as left-to-right (right-to-left). We note 
that the condition for the transformation describing the mapping between S and T matrices to 
be non- singular is 

S13S24 - s 23 si4 7^ 0. (S25) 

This condition is violated when the couplings Mi, M 2 (Fig. SI) are in phase and hence a 
transfer matrix cannot be defined in such a case. 

Similarly, the scattering matrix of the phase shifting stage, 
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yields a transfer matrix of the form 

Tps - 
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Finally, by exploiting the fact that down-conversion is just the inverse operation of the phe- 
nomenon of up-conversion we obtain the transfer matrix of the second UDC stage as 



T DC = Fx T„£. x F, 



(S28) 



where F is the flip matrix of the form F = ia x <8> h required to preserve the consistency in 
labelling the input and output amplitudes. The forward and backward propagation diagrams 
in Fig. 3b of the main text have been calculated using the matrices T uc [Eq. (S24)] and T DC 
[Eq. (S28)]. Each arm of the propagation diagram notes the net forward going amplitude at 
that frequency: for instance, the amplitudes for first UDC stage performing upconversion in 
the forward propagation diagram (green) can be evaluated using the first and third rows of the 
matrix in Eq. (S24) with = 0. The net forward going amplitude contributed by a' to ct+ 
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can be calculated as t n — t 12 = (t +: LR — t* + RL ) = — i/2 3 / 4 while that due to Oq is given by 
£i3 — tiA = {—tt+,LR — H*+rl) = — 1/2 3 ^ 4 , evaluated for the optimal parameter values 5 = 0, 
5± = l/\/2, a = 1/2 3 / 4 reported in the main text. For the backward propagation, a similar 
exercise using the second and fourth rows of T uc (and T DC ) leads to the propagation diagram 
shown in the lower panel of Fig. 3b in red. 
The total transfer matrix 



Ttotai = T DC2 {<P = tt/4) x T PS (9 = it/ 2) x T UCl (<l> = 0) 



(S29) 



establishes the relationship between the low frequency signals on the left and the right ports of 
the cascade as 
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The matrix in Eq. (S30) has been evaluated for the resonant case 8 = and the parameter 
values 

6 ± = a = 2^4- ( S32 ) 

On evaluating the total scattering matrix Stotai from T total , we obtain the four port circulator 
reported in Eq. (9) of the main text. The device behaves nonreciprocally since S T ^ S [1]. 



Modulation ellipse. This scheme, providing a geometrical visualization of the action of our 
device, represents the superposition of two sinusoidal signals as an ellipse in the plane defined 

— * 

by the quadratures / and Q. In general for any two complex phasors a and b rotating in opposite 
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directions 

ae^ + be~ ilJt = Re[(a + b*)e wt ] + ilm[(a - b*)e iu}t ] . (S33) 

S v ' S - ' 

/ Q 

The magnitudes and phases of the two complex signals (4 quantities in total) are encoded as 
different properties of a colored ellipse in the IQ plane: the semi-major axis of the ellipse 
equals p + = \a\ + \b\ while the semi minor axis equals p_ = | \a\ — \b\ |, the angle with the / 
axis equals (9 a — 9b)/2 and the location of the colors on the ellipse represents the phase angle 

(0 a + 9 b )/2. 

The output at each stage in the active circulator comprises two modes - the spatially orthog- 
onal modes (a' , a^) and the sidebands (a+, a_) that can be represented as phasors with opposite 
sense of rotations in a frame rotating at the carrier frequency u c . Thus we can faithfully map the 
detailed dynamics of the device after each stage using the modulation ellipse representation and 
plot the combined output signal in the IQ plane with / = Re[(a+6)e lw *] and Q = Im[(a— b)e wt ] 
(a and b representing the relevant modes at each stage). 

Now we present examples of two different kinds of phase rotations and the resultant trans- 
formations ("rotations") of the modulation ellipse (see Fig. S3 for additional examples on 
modulation ellipse representation). 

1 . Phase shift: 

The action of a phase shifter which performs frequency independent phase rotations of 
both the phasors can be described using the transformations: 

a^ae 10 and b^be id . (S34) 

On using the above and performing the analysis in the IQ plane, we obtain the expressions 
for new coordinates as 

/ = Re[ae l(ult+e) + Ve^ 1 ^} (S35) 
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Q = lm[ae iiu)t+9) - b*e l{u)t -% (S36) 

The action of such an operation can be easily visualized using the modulation ellipse as 
shown in Fig. S2(e). 

2. Free evolution: 

In contrast to the transformation described above, we now consider the rotation preformed 
by a mere time evolution of the two counter-rotating phasors (say by passage through a 
transmission line). In such a case, the phases of the two signals continue to evolve in 
opposite directions collecting a phase 5 in time t (5 = cut), 

a^ae +iS and b ^ be~ iS . (S37) 

The IQ coordinates are calculated as 

I = Re[{a + b*)e t(u}t+5) ] (S38) 
Q = lm[(a-b*)e i ^ t+5) ]. (S39) 

It is immediately evident, that under time evolution, there is only a trivial phase accumu- 
lation leading to change of relative positions of the two phasors along the circumference 
of the modulation ellipse with no rigid rotation of the ellipse, Fig. S2(f). 

In Fig. 3C of the main text, we have used modulation ellipses to represent the dynamics at 
each stage of the active circulator. In the forward propagation direction (L to R), the output 
ellipse representing ports 3 and 4 has the same orientation as the input ellipse representing ports 
1 and 2; only the average phase of the total output changes as indicated by the change of color 
along the circumference of the output ellipse. In the backward propagation direction (R to L), 
however, the ellipse obtained at ports 1 and 2 is rotated by ir/2 with respect to the input ellipse 
of ports 3 and 4, in addition to the trivial average phase change. This indicates the swapping of 
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Figure S2: Modulation ellipse representation of two phasors. In each of the panels, the first 
column describes the phasors under consideration, the second column gives a precise mathemat- 
ical formula for them and the third column shows the corresponding modulation ellipse. In a-d, 
we show both the input modes and the resultant modulation ellipse. The tilt of the ellipse with 
respect to the / axis represents the relative phase between the two modes [(9 a — 9b)/ 2] while the 
color along the ellipse represents the average initial phase of the two modes [(9 a + 9 b )/2], with 
yellow representing the position of (or 2%). Figs, (e) and (f) represent the resulting ellipses 
on performing the indicated transformations on the ellipse in (d). 
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the transmitted amplitudes - port 1 (port 2) receiving the input at port 4 (port 3) - that leads to 
nonreciprocal transmission characteristics of the device. 
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